The phylogenomic landscape of extended-spectrum ß-lactamase producing Citrobacter species isolated from surface water.

BACKGROUND: Citrobacter species are Gram-negative opportunistic pathogens commonly reported in nosocomial-acquired infections. This study characterised four Citrobacter species that were isolated from surface water in the North West Province, South Africa. RESULTS: Phenotypic antimicrobial susceptibility profiles of the isolates demonstrated their ability to produce the extended-spectrum ß-lactamase (ESBL). Whole genomes were sequenced to profile antibiotic resistance and virulence genes, as well as mobile genetic elements. In silico taxonomic identification was conducted by using multi-locus sequence typing and average nucleotide identity. A pangenome was used to determine the phylogenomic landscape of the Citrobacter species by using 109 publicly available genomes. The strains S21 and S23 were identified as C. braakii, while strains S24 and S25 were C. murliniae and C. portucalensis, respectively. Comparative genomics and sequenced genomes of the ESBL-producing isolates consisted of n = 91; 83% Citrobacter species in which bla-CMY-101 (n = 19; 32,2%) and bla-CMY-59 (n = 12; 38,7%) were prevalent in C. braakii, and C. portucalensis strains, respectively. Macrolide (acrAB-TolC, and mdtG) and aminoglycoside (acrD) efflux pumps genes were identified in the four sequenced Citrobacter spp. isolates. The quinolone resistance gene, qnrB13, was exclusive to the C. portucalensis S25 strain. In silico analysis detected plasmid replicon types IncHI1A, IncP, and Col(VCM04) in C. murliniae S24 and C. portucalensis S25, respectively. These potentially facilitate the T4SS secretion system in Citrobacter species. In this study, the C. braakii genomes could be distinguished from C. murliniae and C. portucalensis on the basis of gene encoding for cell surface localisation of the CPS (vexC) and identification of genes involved in capsule polymer synthesis (tviB and tviE). A cluster for the salmochelin siderophore system (iro-BCDEN) was found in C. murliniae S24. This is important when it comes to the pathogenicity pathway that confers an advantage in colonisation. CONCLUSIONS: The emerging and genomic landscapes of these ESBL-producing Citrobacter species are of significant concern due to their dissemination potential in freshwater systems. The presence of these ESBL and multidrug-resistant (MDR) pathogens in aquatic environments is of One Health importance, since they potentially impact the clinical domain, that is, in terms of human health and the agricultural domain, that is, in terms of animal health and food production as well as the environmental domain.

Genetic diversity and whole genome sequence analysis data of multidrug resistant atypical enteropathogenic Escherichia coli O177 strains: An assessment of food safety and public health implications.

Atypical enteropathogenic E. coli (aEPEC) strains are emerging pathogens responsible for fatal diarrhoea in humans worldwide. The purpose of this study was to investigate genetic diversity, virulence and antimicrobial resistance profiles of aEPEC O177 strains isolated from faeces of cattle reared in intensive and extensive production systems in South Africa. A total of 96 multidrug resistant (MDR) aEPEC O177 isolates were typed using enterobacterial repetitive intergenic consensus (ERIC) and random amplified polymorphism DNA (RAPD) typing. The resistome, virulome and mobilome of two aEPEC O177 isolates were investigated using WGS analysis. The ERIC typing was efficient and reproducible with a discriminatory index of 0.95. RAPD typing had poor reproducibility with satisfactory discriminatory power of 0.859. The dendrograms constructed based on ERIC and RAPD banding patterns produced 9 and 8 clusters, respectively, which indicate genetic variation among E. coli O177 isolates. WGS analysis revealed that CF-154-A and CF-335-B) isolates belonged to the O177 serotype with H7 and H21, respectively. Both isolates harboured several virulome genes such as intimin (eaeA), haemolysin (hlyA and hlyE), translocated iron receptor (tir), Type III secretion system (eprH, gspL and prgH), bssR and bssS. However, genes encoding shiga toxins were not found in either isolate. Antibiotic resistance genes such as ampC, tet, ermB, sul2, strB AcrD, aph(6)-Ic, aph(6)-Ib, aph(3″)-I, ant (3″)-1a AcrA and acrE were found in the E. coli O177 strains. Furthermore, genome annotation results indicated that both isolates carried plasmids, insertion sequences, prophages and cluster of regularly interspaced short palindromic repeats (CRISPR) type I. Based on in silico multi locus typing (MLST) analysis, the two isolates were assigned to different sequence types (CF-154-A, ST-1308 and CF-335-B, ST-58). Whole genome multi locus typing tree showed that our isolates clustered with E. coli O177:H21 (reference), suggesting the close genomic relatedness among the strains. Overall, these findings showed that cattle carry genetically diverse E. coli O177 strains, which harbour a repertoire of virulome, resistome and mobilome genes. This highlights a need for multidrug resistant E. coli O177 strain surveillance in cattle.

Screening of Rhizosphere Bacteria and Nematode Populations Associated with Soybean Roots in the Mpumalanga Highveld of South Africa.

Soybean is among South Africa's top crops in terms of production figures. Over the past few years there has been increasingly more damage caused to local soybean by plant-parasitic nematode infections. The presence of Meloidogyne (root-knot nematodes) and Pratylenchus spp. (root lesion nematodes) in soybean fields can cripple the country's production, however, little is known about the soil microbial communities associated with soybean in relation to different levels of Meloidogyne and Pratylenchus infestations, as well as the interaction(s) between them. Therefore, this study aimed to identify the nematode population assemblages and endemic rhizosphere bacteria associated with soybean using Next Generation Sequencing (NGS). The abundance of bacterial genera that were then identified as being significant using linear discriminant analysis (LDA) Effect Size (LEfSe) was compared to the abundance of the most prevalent plant-parasitic nematode genera found across all sampled sites, viz. Meloidogyne and Pratylenchus. While several bacterial genera were identified as significant using LEfSe, only two with increased abundance were associated with decreased abundance of Meloidogyne and Pratylenchus. However, six bacterial genera were associated with decreased Pratylenchus abundance. It is therefore possible that endemic bacterial strains can serve as an alternative method for reducing densities of plant-parasitic nematode genera and in this way reduce the damages caused to this economically important crop.

Wastewater treatment works change the intestinal microbiomes of insectivorous bats.

Mammals, born with a near-sterile intestinal tract, are inoculated with their mothers' microbiome during birth. Thereafter, extrinsic and intrinsic factors shape their intestinal microbe assemblage. Wastewater treatment works (WWTW), sites synonymous with pollutants and pathogens, receive influent from domestic, agricultural and industrial sources. The high nutrient content of wastewater supports abundant populations of chironomid midges (Diptera), which transfer these toxicants and potential pathogens to their predators, such as the banana bat Neoromicia nana (Vespertilionidae), thereby influencing their intestinal microbial assemblages. We used next generation sequencing and 16S rRNA gene profiling to identify and compare intestinal bacteria of N. nana at two reference sites and two WWTW sites. We describe the shared intestinal microbiome of the insectivorous bat, N. nana, consisting of seven phyla and eleven classes. Further, multivariate analyses revealed that location was the most significant driver (sex, body size and condition were not significant) of intestinal microbiome diversity. Bats at WWTW sites exhibited greater intestinal microbiota diversity than those at reference sites, likely due to wastewater exposure, stress and/or altered diet. Changes in their intestinal microbiota assemblages may allow these bats to cope with concomitant stressors.

Antifungal agents, yeast abundance and diversity in surface water: Potential risks to water users.

South African surface waters are subject to various forms of pollution. Recent findings in aquatic systems suggest an association exists between yeast diversity, chemical pollutants and land coverage, which are important water quality determinants. Yeast abundance and diversity, as well as antifungal agents in two river systems in South Africa, were investigated and related to the existing land coverage. Yeast abundance and diversity were determined from environmental DNA by quantitative polymerase chain reaction and next-generation sequencing, respectively, of the 26S ribosomal ribonucleic acid (rRNA) gene. Antifungal agents were qualitatively and/or quantitatively detected by ultra-high-pressure liquid chromatography-mass spectrometry. Analyses of 2 031 714 high-quality 26S rRNA sequences yielded 5554 amplicon sequence variants (ASVs)/species. ASV richness and Shannon-Wiener index of diversity reflected the southward flow of the river with higher values observed downstream compared to the upstream. Fluconazole concentrations were quantifiable in only two samples; 178 and 271 ng L-1. Taxonomically, at least 20 yeast species were detected, including the dominant Candida tropicalis, Cryptococcus spp. as well as the lesser dominant Bensingtonia bomiensis, Fereydounia khargensis, Hericium erinaceus, Kondoa changbaiensi, Pseudozyma spp. and Sphacelotheca pamparum. The two dominant species are known opportunistic pathogens which had antifungal resistant traits in previous studies from the same rivers and therefore is a public health threat. The present study provides further evidence that yeasts should be included as part of water quality parameters, especially in developing countries where much of the population are economically disadvantaged, and also immunocompromised due to age and disease.

Complete Genomic Analysis of VRE From a Cattle Feedlot: Focus on 2 Antibiotic Resistance.

Practices in intensive animal farming such as the extensive use of antimicrobials have significant impacts on the genetic make-up of bacterial communities, especially on that of human/animal commensals. In this report, whole genome sequencing of two vancomycin-resistant enterococci (VRE) isolates from a cattle feedlot in the North West Province, South Africa, was used to highlight the threats that extensive antimicrobial usage in intensive animal rearing represents for environmental microbiomes and the food chain. The genomic DNA of the studied strains was extracted using a DNA extraction kit. Whole-genome sequencing was performed through next-generation sequencing. The genomes of Enterococcus durans strain NWUTAL1 and Enterococcus gallinarum strain S52016 consisted of 3,279,618 and 2,374,946 bp, respectively with G + C contents of 40.76 and 43.13%, respectively. Antibiotic resistance genes (ARG), plasmids and virulence factors (involved in biofilm formation, colonization and copper/silver efflux system), were detected in the genomes of both strains. The presence of these genetic determinants in the studied strains is a cause for concern as they may disseminate and find their way into the food chain via horizontal gene transfer amongst bacteria of the different ecological niches. Issues of this nature cannot be undermined and are relevant as far as food safety is concerned.

Whole-genome sequence of multi-drug resistant Pseudomonas aeruginosa strains UY1PSABAL and UY1PSABAL2 isolated from human broncho-alveolar lavage, Yaoundé, Cameroon.

Pseudomonas aeruginosa has been implicated in a wide range of post-operation wound and lung infections. A wide range of acquired resistance and virulence markers indicate surviving strategy of P. aeruginosa. Complete-genome analysis has been identified as efficient approach towards understanding the pathogenicity of this organism. This study was designed to sequence the entire genome of P. aeruginosa UY1PSABAL and UY1PSABAL2; determine drug-resistance profiles and virulence factors of the isolates; assess factors that contribute toward stability of the genomes; and thereafter determine evolutionary relationships between the strains and other isolates from similar sources. The genomes of the MDR P. aeruginosa UY1PSABAL and UY1PSABAL2 were sequenced on the Illumina Miseq platform. The raw sequenced reads were assessed for quality using FastQC v.0.11.5 and filtered for low quality reads and adapter regions using Trimmomatic v.0.36. The de novo genome assembly was made with SPAdes v.3.13 and annotated using Prokka v.2.1.1 annotation pipeline; Rapid Annotation using Subsytems Technology (RAST) server v.2.0; and PATRIC annotation tool v.3.6.2. Antimicrobial resistance genes and virulence determinants were searched through the functional annotation data generated from Prokka, RAST and PATRIC annotation pipelines; In addition to ResFinder and Comprehensive Antibiotic Resistance Database (CARD) which were employed to determine resistance genes. The PHAge Search Tool Enhanced Release (PHASTER) web server was used for the rapid identification and annotation of prophage sequences within bacterial genome. Predictive secondary metabolites were identified with AntiSMASH v.5.0. Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR) and cas genes regions were also investigated with the CRISPRone and CRISPRFinder server. The genome sizes of 7.0 and 6.4 Mb were determined for UY1PSABAL and UY1PSABAL2 strains with G+C contents of 66.1% and 66.48% respectively. ß-lactamines resistance genes blaPAO, aminoglycoside phosphorylating enzymes genes aph(3')-IIb, fosfomycine resistance gene fosA, vancomycin vanW and tetracycline tetA were among identified resistance genes harboured in both isolates. UY1PSABAL bore additional aph(6)-Id, aph(3'')-Ib, ciprofloxacin-modifying enzyme crpP and ribosomal methylation enzyme rmtB. Both isolates were found harbouring virulence markers such as flagella and type IV pili; and also present various type III secretion systems such as exoA, exoS, exoU, exoT. Secondary metabolites such as pyochelin and pyoverdine with iron uptake activity were found within the genomes as well as quorum-sensing systems, and various fragments for prophages and insertion sequences. Only the UY1PSABAL2 contains CRISPR-Cas system. The phylogeny revealed a very close evolutionary relationship between UY1PSABAL and the similar strain isolated from Malaysia; the same trend was observed between UY1PSABAL2 and the strain from Chinese origin. Complete analyses of the entire genomes provide a wide range of information towards understanding pathogenicity of the pathogens in question.

Combining physicochemical properties and microbiome data to evaluate the water quality of South African drinking water production plants.

Anthropogenic activities in catchments used for drinking water production largely contaminates source waters, and this may impact the quality of the final drinking water product. These contaminants may also affect taxonomic and functional profiles of the bacterial communities in the drinking water. Here, we report an integrated insight into the microbiome and water quality of four water treatment plants (NWC, NWE, WCA and NWG) that supply portable water to communities in South Africa. A new scoring system based on combined significant changes of physicochemical parameters and microbial abundance from raw to treated water was used to evaluate the effectiveness of the treatment plants at water purification. Physicochemical parameters which include total soluble solids, turbidity, pH, nitrites and phosphorus among others, were measured in source, treated, and distributed water. There were general statistically significant (P ≤ 0.05) differences between raw and treated water, demonstrating the effectiveness of the purification process. Illumina sequencing of the 16S rRNA gene was used for taxonomic profiling of the microbial communities and this data was used to infer functional attributes of the communities. Structure and composition of the bacterial communities differed significantly (P < 0.05) among the treatment plants, only NWE and NWG showed no significant differences (P > 0.05), this correlated with the predicted functional profile of the microbial communities obtained from Phylogenetic Investigation of Communities by Reconstruction of Observed States (PICRUSt), as well as the likely pollutants of source water. Bacteroidetes, Chlorobi and Fibrobacteres significantly differed (P < 0.05) between raw and distributed water. PICRUSt inferred a number of pathways involved in the degradation of xenobiotics such as Dichlorodiphenyltrichloroethane, atrazine and polycyclic aromatic hydrocarbons. More worryingly, was the presence of pathways involved in beta-lactam resistance, potential pathogenic Escherichia coli infection, Vibrio cholerae infection, and Shigellosis. Also present in drinking and treated water were OTUs associated with a number of opportunistic pathogens.

Inside environmental Clostridium perfringens genomes: antibiotic resistance genes, virulence factors and genomic features.

Until recently, research has focused on Clostridium perfringens in clinical settings without considering environmental isolates. In this study, environmental genomes were used to investigate possible antibiotic resistance and the presence of virulence traits in C. perfringens strains from raw surface water. In silico assembly of three C. perfringens strains, DNA generated almost complete genomes setting their length ranging from 3.4 to 3.6 Mbp with GC content of 28.18%. An average of 3,175 open reading frames was identified, with the majority associated with carbohydrate and protein metabolisms. The genomes harboured several antibiotic resistance genes for glycopeptides, macrolide-lincosamide-streptogramin B, ß-lactam, trimethoprim, tetracycline and aminoglycosides and also the presence of several genes encoding for polypeptides and multidrug resistance efflux pumps and 35 virulence genes. Some of these encode for haemolysins, sialidase, hyaluronidase, collagenase, perfringolysin O and phospholipase C. All three genomes contained sequences indicating phage, antibiotic resistance and pathogenic islands integration sites. A genomic comparison of these three strains confirmed high similarity and shared core genes with clinical C. perfringens strains, highlighting their health security risks. This study provides a genomic insight into the potential pathogenicity of C. perfringens present in the environment and emphasises the importance of monitoring this niche in the future.

Draft Genome Sequences of Potentially Pathogenic Clostridium perfringens Strains from Environmental Surface Water in the North West Province of South Africa.

Surface water systems in South Africa are experiencing a major decline in quality due to various anthropogenic factors. This poses a possible health risk for humans. Here, we present the draft genome sequences of three Clostridium perfringens isolates obtained from a fecally polluted river system in the North West province of South Africa.

Draft Genome Assemblages of 10 Xanthomonas vasicola pv. zeae Strains, Pathogens Causing Leaf Streak Disease of Maize in South Africa.

Maize bacterial leaf streak disease has spread across maize crops in South Africa and therefore potentially poses a threat to maize production and food security. Until recently, this pathogen was identified as a Xanthomonas campestris pathovar, whereas our South African genomes seem to be more divergent and create their own subclade.

Ecological guild and enzyme activities of rhizosphere soil microbial communities associated with Bt-maize cultivation under field conditions in North West Province of South Africa.

Insecticidal proteins expressed by genetically modified Bt maize may alter the enzymatic and microbial communities associated with rhizosphere soil. This study investigated the structure and enzymatic activity of rhizosphere soil microbial communities associated with field grown Bt and non-Bt maize. Rhizosphere soil samples were collected from Bt and non-Bt fields under dryland and irrigated conditions. Samples were subjected to chemical tests, enzyme analyses, and next generation sequencing. Results showed that nitrate and phosphorus concentrations were significantly higher in non-Bt maize dryland soils, while organic carbon was significantly higher in non-Bt maize irrigated field soil. Acid phosphatase and ß-glucosidase activities were significantly reduced in soils under Bt maize cultivation. The species diversity differed between fields and Bt and non-Bt maize soils. Results revealed that Actinobacteria, Proteobacteria, and Acidobacteria were the dominant phyla present in these soils. Redundancy analyses indicated that some chemical properties and enzyme activities could explain differences in bacterial community structures. Variances existed in microbial community structures between Bt and non-Bt maize fields. There were also differences between the chemical and biochemical properties of rhizosphere soils under Bt and non-Bt maize cultivation. These differences could be related to agricultural practices and cultivar type.

Culicoides species composition and environmental factors influencing African horse sickness distribution at three sites in Namibia.

African horse sickness (AHS) is one of the most lethal infectious, non-contagious, vector-borne disease of equids. The causative agent, African horse sickness virus (AHSV) is transmitted via Culicoides midges (Diptera: Ceratopogonidae). AHS is endemic to Namibia but detailed studies of Culicoides communities and influencing environmental parameters are limited. This study aims to determine the Culicoides species composition at three different sites and to assess environmental parameters influencing the geographical distribution of AHS in Namibia. Weekly collections of Culicoides were made during the AHS peak season from January to May for 2013 and 2014 using the Onderstepoort 220V UV-light trap. Out of 397 collections made, 124 collections (3287 Culicoides) were analysed for AHSV presence with RT-qPCR. A total of 295 collections were analysed for total Culicoides (all collected Culicoides individuals) and in 75% of these collections the Culicoides were identified to species level. C. imicola was the dominant species with proportional representation of 29.9%. C. subschultzei, C. exspectator and C. ravus each contribute more than 10% to the species composition. The lowest number of Culicoides was collected at Aus 9980, a total of 21819 at Windhoek and the highest number at Okahandja 47343. AHSV was present at all three sites during 2013 but only in Windhoek and Okahandja during 2014. Multivariate analyses of data from the two year survey indicate the environmental parameters in order of importance for the distribution of AHS in Namibia as precipitation>temperature>clay>relative humidity>NDVI. The implication of these findings is that any precipitation event increases Culicoides numbers significantly. Together with these results the high number of species found of which little is known regarding their vector competence, add to the complexity of the distribution of AHS in Namibia.

Detection of African horse sickness virus in Culicoides imicola pools using RT-qPCR.

African horse sickness (AHS) is an infectious, non-contagious arthropod-borne disease of equids, caused by the African horse sickness virus (AHSV), an orbivirus of the Reoviridae family. It is endemic in sub-Saharan Africa and thought to be the most lethal viral disease of horses. This study focused on detection of AHSV in Culicoides imicola (Diptera: Ceratopogonidae) pools by the application of a RT-qPCR. Midges were fed on AHSV-infected blood. A single blood-engorged female was allocated to pools of unfed nulliparous female midges. Pool sizes varied from 1 to 200. RNA was extracted and prepared for RT-qPCR. The virus was successfully detected and the optimal pool size for the limit of detection of the virus was determined at a range between 1 to 25. Results from this investigation highlight the need for a standardized protocol for AHSV investigation in Culicoides midges especially for comparison among different studies and for the determination of infection rate.

PCR-denaturing gradient gel electrophoresis analysis of microbial community in soy-daddawa, a Nigerian fermented soybean (Glycine max (L.) Merr.) condiment.

Soy-daddawa, a fermented soybean (Glycine max (L.) Merr.) condiment, plays a significant role in the culinary practice of West Africa. It is essential to understand the microbial community of soy-daddawa for a successful starter culture application. This study investigated the microbial community structure of soy-daddawa samples collected from Nigerian markets, by PCR-denaturing gradient gel electrophoresis (DGGE) targeting the V3-V5 region of the 16S rRNA gene of bacteria and internal transcribed spacer 2 (ITS2) region of fungi. Six bacterial and 16 fungal (nine yeasts and seven molds) operational taxonomic units (OTUs)/species were obtained at 97% sequence similarity. Taxonomic assignments revealed that bacterial OTUs belonged to the phyla Firmicutes and Actinobacteria, and included species from the genera Atopostipes, Bacillus, Brevibacterium and Nosocomiicoccus. Densitometric analysis of DGGE image/bands revealed that Bacillus spp. were the dominant OTU/species in terms of population numbers. Fungal OTUs belonged to the phyla Ascomycota and Zygomycota, and included species from the genera, Alternaria, Aspergillus, Candida, Cladosporium, Dokmaia, Issatchenkia, Kodamaea, Lecythophora, Phoma, Pichia, Rhizopus, Saccharomyces and Starmerella. The majority of fungal species have not been previously reported in soy-daddawa. Potential opportunistic human pathogens such as Atopostipes suicloacalis, Candida rugosa, Candida tropicalis, and Kodamaea ohmeri were detected. Variation in soy-daddawa microbial communities amongst samples and presence of potential opportunistic pathogens emphasises the need for starter culture employment and good handling practices in soy-daddawa processing.

Evolutionary history of Phakopsora pachyrhizi (the Asian soybean rust) in Brazil based on nucleotide sequences of the internal transcribed spacer region of the nuclear ribosomal DNA

Phakopsora pachyrhizi has dispersed globally and brought severe economic losses to soybean growers. The fungus has been established in Brazil since 2002 and is found nationwide. To gather information on the temporal and spatial patterns of genetic variation in P. pachyrhizi, we sequenced the nuclear internal transcribed spacer regions (ITS1 and ITS2). Total genomic DNA was extracted using either lyophilized urediniospores or lesions removed from infected leaves sampled from 26 soybean fields in Brazil and one field in South Africa. Cloning prior to sequencing was necessary because direct sequencing of PCR amplicons gave partially unreadable electrophoretograms with peak displacements suggestive of multiple sequences with length polymorphism. Sequences were determined from four clones per field. ITS sequences from African or Asian isolates available from the GenBank were included in the analyses. Independent sequence alignments of the ITS1 and ITS2 datasets identified 27 and 19 ribotypes, respectively. Molecular phylogeographic analyses revealed that ribotypes of widespread distribution in Brazil displayed characteristics of ancestrality and were shared with Africa and Asia, while ribotypes of rare occurrence in Brazil were indigenous. The results suggest P. pachyrhizi found in Brazil as originating from multiple, independent long-distance dispersal events.